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1. Introduction 



The first part of these proceedings will give an update on our results for the lowest moments of 
the leading twist meson distribution amplitudes (DAs) |[l], §]. Distribution amplitudes contain non- 
perturbative QCD effects that appear e.g. in hard exclusive processes but are universal hadronic 
properties and so do not depend on the process itself. They are important for form factors at large 
q 2 or B-decays and can be related to the Bethe-Salpeter wave function. Here we provide results 
for the pseudo-scalar DAs of K and % as well as for the vector DAs of K*, p and (p. More details 
on these calculations by the RBC and UKQCD collaborations will be presented in a forthcoming 
paper. 

The second part is devoted to the non-perturbative renormalisation of quark bilinears with and 
without derivatives. Renormalisation is necessary to obtain physical results for e.g. the matrix ele- 
ments of moments of DAs or decay constants. We follow the Rome-Southampton method [^] with 
momentum sources [Q] to calculate the renormalisation constants. The first is aimed at reducing 
uncertainties in the perturbative renormalisation while the latter is to improve upon the statistical 
errors from the standard point source approach. 

The calculations are done on lattices with 24 3 x 64 x 16 and 16 3 x 32 x 16 points at a lattice 
spacing of a -1 = 1.729(28) GeV. The ensembles have been generated by the RBC/UKQCD col- 
laborations using Nf = 2 + 1 domain wall fermions with an Iwasaki gauge action. Our light quark 
masses range from am q = 0.005 to 0.03 corresponding to pion masses from 331MeV to 672 MeV. 
The strange quark mass is kept fixed at am s = 0.04. Details on the ensembles have been reported 
in©!]. 

2. Meson distribution amplitudes 

The meson DAs are defined as non-local matrix elements on the light-cone. The leading twist 
pseudo-scalar and (longitudinal) vector DAs are 

(0|^ Fl (z)y M y 5 ^(z,-z)y^(-z)^^ (2.1) 
<0| y Fl (z)y^(z,-z)Y Fl (-z) \VM)\^=ftm v pfi^ j\e { ^4^^). (2.2) 

Here B, = 2x — I with x and 1 — x the momentum fractions of the two quarks. The Wilson line 
% ensures gauge invariance and the quark flavours F, are chosen to match the pseudo-scalar and 
vector mesons. All of the DAs 0(^,ju) are normalised to unity when integrated over £ G [—1,1]. 
On the lattice we only access moments of the DAs, 

(D(M) = /d^'X^At), (2.3) 

which appear in matrix elements of local operators with n derivatives, see e.g. [Q]. We extract the 
bare values for 312 (n = 1,2) from the lattice using ratios of two-point functions. To give one 
example, let us consider the second moment of the (longitudinal) vector meson DA, /(£ 2 \"' bare Qj), 
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For large Euclidean times t and T — t, the conelation functions give for a ratio like 



i<^ 2 >"' bai ' e tanh((f-r/2)£ v 





^Wv}(*,OkJ(o)|o} 







1 ( , ^PpPjiPvPo 



X - -gpaPuPv ~ g^aPpPv ~ gvaPpP^i H 5 " ( 2 " 4 ) 

-5 V mr, 



Here Z?y and my are the energy and mass of the vector meson with interpolating field V. The 
operator with n derivatives is given by 

^Wi-M(*>0 = ¥F l (x,t)Y{nD fll ...D fln} YF 2 (x,t), (2.5) 



where the {. . . } denote symmetrisation of the indices and subtraction of traces. Note that by choos- 
ing directions such as e.g. pi = G = 2, v = 4 and p = 1 one unit of momentum p\ 7^ is enough to 
extract (^IIMe ^ from Eq q 

The results of fits to ratios like Eq. (2A) are plotted for both lattice sizes in Fig. RL Also shown 
are linear extrapolations of these bare results to the chiral limit along with their error bands. For 
the lowest moment of the kaon DA such a linear extrapolation in (m s — m q ) is predicted by NLO 
chiral perturbation theory [JSj] . Note that we have to account for a slightly wrong strange quark 
mass of am & = 0.04 in our simulation instead of the physical am s = 0.0343(16) [EJ. In case of the 
lowest moment for K and K*, this is done [§] by extrapolating in (ra s — m q ) where we can simply 
enter the correct physical quark masses. For the remaining second moments, a correction in m s 
could be obtained by comparing the results for mesons with and without strange quarks. However, 
within errors these differences are small enough to neglect any corrections. Fig. [I] also shows that 
finite size effects are small and not significant within errors. Finally, the preliminary results for our 
((§") (n = 1,2) at the physical point are given in Tab. [j], along with results by other groups where 
available. Our values are quoted using perturbative renormalisation to the MS scheme at a scale 
of pi = 2GeV. We have already pointed out the clearly observable SU (3)-breaking effects for the 
kaon DA in [jl]]. The same now holds for the K*. 



Table 1: Our preliminary results at the physical point and both our lattice sizes, in MS at jl — 2GeV using 
perturbative renormalisation. Also included are results from the given references (at jl = 2GeV). 











ml 


24 3 


0.02893(87)(166) 


0.267(11)(13) 


0.0342(16)(21) 


0.248(17)(12) 


16 3 


0.0277(17)(16) 


0.282(17)(14) 


0.0297(11)(16) 


0.255(13)(13) 




0.0272(5) 


0.260(6) 


0.033(2)(4) 














24 3 


0.272(15)(13) 


0.237(36)(12) 


0.246(10)(12) 




16 3 


0.274(34)(13) 


0.245(27)(12) 


0.245(11)(12) 




[§ 0.269(39) 
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Figure 1: Linear extrapolations of the various mo- 
ments of the DAs for pseudo-scalar and vector mesons 
against the mass difference m s — m q or the quark mass. 
Shown are the bare results for K, K*, K, p and (j) for 
both lattice sizes. The solid vertical lines mark the 
physical point with the uncertainty in m s in case of the 
lowest moments of K and K*. 



3. Non-perturbative renormalisation 

The previous section made use of perturbative renormalisation only. Since this introduces 
uncertainties due to the known bad convergence of the perturbative expansion, our task now is to 
compute the necessary renormalisation constants non-perturbatively. For this, we use the Rome- 
Southampton method which employs a simple renormalisation condition that is useful for any 
regularisation [gj: 



K ff {p)=Z e ^)Z- l h^{p) 



1. (3.1) 
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Here A^ are \p) is the renormalised (bare) vertex amplitude (definition in Eq. (|3.3[)), Z q is the 
quark field renormalisation (y = Z^ 2 \j/ hme ) and Z@(\i) the desired renormalisation constant (& = 
Z^^ bare ) at the scale \i. The renormalisation condition should be applied for scales within an 



appropriate window [gg, Aqcd <C \l <C I /a, that we know from our previous calculation Q10|]. 

To obtain the bare vertex function, we start by calculating the unamputated Green's function 
of the operator G between external off-shell quarks with exceptional momenta, i.e. p = p' 

Gffip) = (y(p)0(PW(p)) =£(y 5 [L y S(x,yyr>f Y 5 J r (x,x') [£ z5 (V,z) e ^]) . (3.2) 

Where the operator is written as 0(q) = Y^xjtf ^p{iqx)\j/(x)Jr(x,x / )Y(x') with Dirac structure T and 
possible derivatives. We have already written the quark propagators and their Fourier transform in 
Eq. ( |3.2[ ) in a way suggesting the use of momentum sources [Q]. Instead of using a point source 
for the inversion, one can perform the inversion on a momentum source e lpz to find a solution for 
S(p)x = HyS(x,y)e ipy ■ We then amputate the Green's function and project onto the bare (tree-level) 
vertex amplitude, 

rW) = (S(p)) G l (G ff (p)) G (S(p)) G l and A,(p) = ^lr{U e {p)P e ) . (3.3) 

Here the subscript G indicates the gauge average and N ensures the overall normalisation of the 
trace. The projector P@ has to match the specific Lorentz and kinematical structure of the operator. 



Examples for quark bilinears without derivatives can be found e.g. in [10|. The operators involving 
derivatives have a more complicated decomposition that is consistent with their symmetries and 
tracelessness. Thus more care has to be taken to project onto the correct tree-level contribution, see 



e.g. Jl lp. One possibility for an operator like YY{n^v}W would be 

Tr [n^ Dv (p)(7n + 7v)] Ip^.vTr [n^ Dv (p)Y P 



/j<V 



(3.4) 



This particular example averages over all possible space-time components of the operator and pro- 
jector. The symmetrisation of the operator is reflected in the first part of the projection, while 
the second part ensures we pick up the part proportional to the tree-level contribution. We use 
p^ = sin(p^) to compensate the kinematic factors. 

Since this is an ongoing project, let us only mention a few important findings concerning the 
advantages of the momentum sources. Fig. ^ compares results from our earlier calculation using 
point sources, see JlO|], with the current results. Of course the two calculations agree, the important 
point is the much smaller statistical error for the momentum sources as shown in the insert of 
Fig. ||. These results have been obtained from one direction of the momentum for each (ap) 2 and 
up to 25 configurations. The point source results on the other hand used 4 separate sources with 75 
configurations each and average multiple directions for one value of (ap) 2 , [1C]. It is thus possible 
to reduce the computational cost while having smaller statistical errors. 

The smaller errors make it possible to see lattice discretisation errors which is shown in Figs. ||- 
|5[ Let us consider quark bilinears without derivatives first. Here, one can imagine additional 
contributions from terms <=c p 2 . Using latticised momenta p again and expanding the sine function, 
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Figure 2: Comparison of the bare vertex ampli- 
tude for the vector and axial current from point and 
momentum sources. Both on 16 3 x 32 lattices for 
am a = 0.03. 
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Figure 3: The bare vertex amplitude for the scalar 
density together with an indicator for discretisation 
errors, 5? (see text). Again 16 3 x 32, am q = 0.03. 
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Figure 4: Normalised bare vertex for the two source 
types, again 16 3 x 32, am q = 0.03. Only momentum 
sources reveal different discretisation errors. 
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Figure 5: The vertex of the one derivative operator 
for two momenta, split up into the different choices 
of Lorentz indices {pv}. Again 16 3 x 32. 



we expect possible discretisation errors to be proportional to y = ^ {inp^/L^ (taking the 
leading correction only). In Fig. ^ we plot the bare vertex amplitude of the scalar density along 
with y. The deviations from the expected smooth behaviour of the vertex is clearly correlated 
with large changes in .y. This becomes even clearer, when we look at a fixed (ap) 2 but momenta 
in different directions. The bare vertex should not depend on the latter. However, the different 
directions result in different values of S fi . Fig. ^ shows again the scalar density, normalised with 
the mean vector and axial vector vertex. Included are all directions to one momentum as used in 
UlQU, sorted according to the corresponding value for S fi . The different discretisation errors are 
only visible for results obtained with momentum sources however. Some examples of these are 
included in the plot (note that their errors are hidden by the symbols). 

Finally, let us mention that the discretisation errors are more difficult for operators with deriva- 
tives. Here one can distinguish different directions of the Lorentz indices of the operator w.r.t. 
the direction of the momentum. This fact is demonstrated for two momenta and the vertex for 
Y7{^D v j\j/ in Fig. ||. In the continuum all directions are equivalent, on the lattice only choices 
which have similar momentum components agree within errors, e.g. the combinations {pv} = 
{23}, {24}. 
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4. Summary 

We have presented an update on our (preliminary) results for distribution amplitudes which 
have been extended to more pseudo-scalar and vector mesons. The SU (3)-breaking effects already 
found for the K have been confirmed for the K*. We did not see a clear sign of finite size effects. 

We have also presented our first findings using momentum sources for a non-perturbative 
renormalisation of our lattice results. Here we see a clear advantage in reducing the statistical 
errors making it possible to better control the effects due to the lattice discretisation. 
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